Statistics of 3-dimensional Lagrangian turbulence 
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We consider a superstatistical model for a Lagrangian tracer particle in a high-Reynolds number 
turbulent flow. The analytical model predictions are in excellent agreement with recent experimen- 
tal data for flow between counter-rotating disks. In particular, the predicted Lagrangian scaling 
exponents Q agree well with the measured exponents reported in [X. Hu et al., PRL 96, 114503 
(2006)]. The model also correctly predicts the shape of acceleration probability densities, correlation 
functions, statistical dependencies between components and explains the fact that enstrophy lags 
behind dissipation. 
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The full understanding of the statistical properties of 
fluid turbulence still remains a challenging problem in 
theoretical physics. In recent years there has been some 
experimental progress in measuring the statistical prop- 
erties of single tracer particles advected by turbulent 
flows P, [3, S, B 0, ll H ■ These experiments^s well 
as direct numerical simulations (DNS) El, [13, [Q, 
have significantly enhanced our knowledge of the stochas- 
tic properties of Lagrangian turbulence (LT). A variety 
of interesting experimental results has recently been pub- 
lished, in particular for the probability densities of accel- 
erations 1 , [ll , velocity differences 3 , correlation func- 
tions [1, 5 1 , conditional expectations H, 0, 0] , the La- 
grangian scaling exponents [1, 0, Q and the corre- 
sponding /(a) spectra obtained by Legendre transforma- 
tion Q. Many of these new experimental data confirm 
early DNS results obtained in |10| . 

The exact theory of LT based on the Navier-Stokes 
equation still very much lags behind the experimental 
progress. Hence it is important to develop simple theo- 
retical models that provide an explanation for the most 
important statistical features of 3-d LT. Previous models 
have been succesfull in explaining e.g. the 1-d measured 
acceleration statistics but fail to explain the recent exper- 
imental data for the 3-d statistics i4] or the new data for 
the Lagrangian scaling exponents @ [§] . It is thus impor- 
tant to develop theoretical models that explain not just 
one but all of the above measured LT phenomena with 
sufficient precision at the same time, using a consistent 
set of parameters. 

In this letter we will introduce such a model, which 
for the first time simultaneously reproduces the mea- 
sured 3-d statistics, correlation functions, statistical de- 
pendencies between components and scaling exponents. 
We will carefully compare its theoretical predictions 
with the available experimental data, obtaining excel- 




lent agreement. Our model is a natural and physi- 
cally plausible extension of previous LT models based on 
superstatistical stochastic^ differential equations (SDEs) 

. Here the term 'superstatis- 
is a superposition of several 
stochastic processes, a fast one as given by the original 
SDE and a slow one for the parameters of the SDE, which 
are regarded as slowly varying random variables describ- 
ing the changing environment of the Lagrangian tracer 
particle. Our model approximates the high-Reynolds 
number limit of a superstatistical extension of the Saw- 
ford model [SilUll, and is more refined than previous 



models [l3.Tl5l 
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17| | by taking into account both a fluc- 
tuating energy dissipation and a fluctuating enstrophy 
surrounding the test particle. 

We obtain predictions for Lagrangian scaling expo- 
nents (and, by Legendre transform, for multifractal 
spectra) that are in very good agreement with the re- 
cent measurements of the Bodenschatz group [8, 9J. The 
agreement seems to be better than for other models, e.g. 
the Lagrangian multifractal turbulence models compared 
with the data in [3]- At the same time our model cor- 
rectly reproduces the measured probability densities of 
single velocity difference and acceleration components as 
well as those of the modulus, it describes correctly the 
fact that the three acceleration components are not sta- 
tistically independent, it gives the correct conditional ac- 
celeration variance, and it explains the fact that the cor- 
relation function for single acceleration components de- 
cays rapidly whereas that of the modulus decays slowly. 
Finally, the model also explains why the fluctuating en- 
strophy lags behind the fluctuating energy dissipation, as 
numerically observed by Yeung and Pope and exper- 
imentally hy Zeff et al. [23|. To the best of our knowl- 
edge, our model is the flrst LT model that simultaneously 
achieves all this. 

To introduce the model, let us denote the velocity of 
a Lagrangian tracer particle embedded in the turbulent 
flow by v(t). We are interested in velocity differences on 
given time scales r, i.e. the quantity 5v{t) = v{t) — v{t + 
r), which is directly measured in various experiments. 
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To obtain a compact notation, in the following we write 
u{t) := Sv{t). Let us consider a linear superstatistical 
SDE for u{t) of the following form: 



On a time scale t satisfying '~f~^ << t << To- the 
probability density of a single component of the tracer 
particle described by eq. ^ is given by the Gaussian 



■= -ru + j:L{t) 



(1) 



Here L{t) is a rapidly fluctuating stochastic process rep- 
resenting force differences in the liquid on a fast time 
scale, and F and S are 3x3 matrices. We approximate 
L(t) by Gaussian white noise. T(t) and S(i) are matrix- 
valued stochastic processes which evolve on a much larger 
time scale than L{t). 

Generally, the class of superstatistical models de- 
scribed by eq. ([T]) is quite large. A first example was 
introduced in [14] and further developed in [H, [O]- Let 
us here consider a more refined model that takes into 
account two very important facts: A fluctuating local 
energy dissipation rate of the environment surrounding 
the test particle and a fluctuating local enstrophy (ro- 
tational energy). We thus consider as a special case of 
eq. (m the local dynamics 



-ju + Bn X u + aL{t) . 



(2) 



We assume that 7 and B are constants, but the noise 
strength a and the unit vector ft describing a temporary 
rotation axis of the particle evolve stochastically on a 
large time scale and T^, respectively. T^. is of the same 
order of magnitude as the integral time scale Tl , whereas 
7^^ is of the same order of magnitude as the Kolmogorov 
time scale t,j. From the above one gets T^j Tl/t,-/ ^ 
R\ >> 1, i.e. the time scale separation between the slow 
and the fast processes, which is a necessary condition 
for the superstatistics approach to work [24|, increases 
proportional to the Taylor scale Reynolds number R\. 
The time scale » describes the average life time 
of a region of given vorticity surrounding the test particle. 

As it is customary in statistical physics, we define a pa- 
rameter /? := 27/(7^, which in equilibrium statistical me- 
chanics corresponds to the inverse temperature, whereas 
in superstatistical turbulence models [3, [H, [iBl it is a 
formal parameter related to a suitable power e" of the 
fluctuating energy dissipation rate e. In the following 
we adopt the choice k — —1, i.e. I3~^ ^ :/^/^(e)^^/^e, 
where v is the kinematic viscosity and (e) the average 
energy dissipation. To further specify our superstatisti- 
cal model we still have to fix the probability density of 
the stochastic process f}{t), which, motivated by the cas- 
cade picture of turbulence and previous successful models 
d, [H, [13, m, H^, , is assumed to be close to a lognor- 
mal distribution 
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(4) 



Note that this result is independent of B and n [27|. In 
the long-term run the variance of this Gaussian will fluc- 
tuate since cr fluctuates. Hence we get a superposition 
of Gaussians with different variance parameter (3~^ , i.e. 
the marginal stationary distribution of our superstatisti- 
cal system is given by 
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(5) 



no matter what the value of B is. This formula, with /(/3) 
being the lognormal distribution, is in excellent agree- 
ment with experimentally measured histograms [J, [3i 0] ■ 
An example is shown in Fig. 1. One obtains good 
fittings of the data in 3] for all time scales r, with 
w ^ (t/tjj)^'^''^. For very small r the acceleration of 
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FIG. 1: Experimentally measured probability distribution 
p{ux) at smallest time scales (data from for r ~ 0.15 ms) 
and a fit by eq. (0) and ^ with = 1.8. 

the particle is given by /t and hence by trans- 

formation of random variables Pa^io-x) — Tpu^(ux), thus 



2tts 



d(3 13-^'^ exp . 



-(lo^ 



P \2 ■ 



2s2 



(6) 

This formula is in good agreement with the results pre- 
sented in P,[l] (s^ = 3.0, see [Tl] for details). 

Let us now check what type of Lagrangian scaling ex- 
ponents Cj" for velocity increments u are predicted by our 
superstatistical model. From eq. (O and eq. ^ one ob- 
tains the moments 



Here m and s are mean and variance parameters. The 
average /3o of the above log-normal distribution is given 
by /3o = m\fw and the variance by = rn^wiw — 1), 
where w := e" . 



(u^) = (j — l)!!r7i 2u;sJ'^. 



(7) 



The notation (j — 1)!! stands for a product of all odd pos- 
itive integers up to j — 1. Assuming simple scaling laws 
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of the form m 



T^, where a and h are so far ar- 



bitrary real numbers, eq. ([7]) imphes 
r''^ . Hence 







(8) 



It is often assumed that the Lagrangian exponent C,2 is 
equal to 1. From C2 = 1 we get a = \h — 1 thus 



(9) 



where, following the notation of Mordant et al. we 
defined a positive parameter hy \^ = —jh. Note that 
we get a formula analogue to Kolmogorov's K62 theory 
[2^ , however the difference is that this formula is directly 
applicable to the Lagrangian dynamics, it needs not to 
be transformed from an Eulerian to a Lagrangian frame 
(as it was done in [2^ for multifractal models) . Our 
formula ([9]) is in very good agreement with the experi- 
mental data presented in [9], see Fig. 2. The agreement 




of the three components of the Lagrangian par- 

ticle. It is given by the superstatistical average 



1 



(10) 

In particular, for small t the distribution of the absolute 
value \d\ = T^^\u\ of acceleration is given by 



p{\d\) = AT:\d\^p{ax,ay,az 
I2 



(11) 



/33/V(/3)e" = ^" d(3. (12) 



The agreement of this formula with experimentally mea- 
sured distributions of the acceleration modulus has been 
checked in [bJ, taking again for /(/3) the lognormal dis- 
tribution. Excellent agreement was found. Note that 
in cq. pH)) the 3-point density is not the product of 1- 
point densities as given by eq. ©, and hence the model 
naturally introduces statistical dependence between the 
acceleration components. 

We may investigate this effect in a quantitative way, 
by studying the ratios i? :— p{aj:, ay)/ {p{ax)p{ay)). From 
our superstatistical 3-d model we obtain the general pre- 
diction 



R = 



r/3^/V(/3)e 



d/3/„°°/3i/2/(/3)e 



yd/3 
(13) 

which is plotted in Fig. 3 for the example of the lognormal 
distribution /(/3). One obtains a diagram that closely 



FIG. 2: Lagrangian scaling exponents Q as measured by Xu et 
al. and as predicted by the superstatistical model for A'^ — 
0.085 (solid line). Other model predictions are also shown 
(from top to bottom): No intermittency at all (dashed line), 
Lagrangian version of She-Leveque model [9l.[2^ (dotted line) 
and multifractal model of Chevillard et al. [y, (dashed 
dotted line). 

is better than for the other theoretical models compared 
with the data in . Whereas all our predicted exponents 
are within the error bars of the experimental data, the 
predictions of the other models are outside the experi- 
mental error bars for j > 5 (dotted and dashed-dotted 
lines in Fig. 2). 

We may also proceed to the multifractal turbulence 
spectra D{h) defined by D{h) = mij{hj + 1 ~ Q) by a 
Legendre transformation. The information contained in 
the D(h) is the same as that in the Q, and hence our pre- 
dicted multifractal spectra obtained by Legendre trans- 
formation are again in good agreement with the experi- 
mental data presented in Q , better than the predictions 
of the other models. 

Next, let us study the multivariate distribution 
p(ux,Uy,Uz) describing the joint probability distribution 
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FIG. 3: The ratio R = p{ax, ay)/{p{ax)p{ay) as predicted 
by the superstatistical model for — 3. i? 7^ 1 indicates 
statistical dependence of acceleration components. 

resembles the experimentally measured data (Fig. 5 in 
0)- 

All kinds of quantities describing the statistical depen- 
dence can be analytically evaluated for our model. For 
example, one obtains the conditional moments (a^lay) as 



= (J-1) 



_J„°°d/3/3^/(/3)e- 



/o°°d/3/3i/2/(/3)e- 



(14) 



4 



for j even (they vanish for j odd). Moreover, for fognor- 
mal superstatistics one obtains 



«4) 
«)(4> 



(15) 



even) which yields a relation between the flatness 
parameter w and the statistical dependencies of the ac- 
celeration components that can be checked in future ex- 
periments. 

Our model also allows for the calculation of temporal 
correlation functions. In particular we may be interested 
in temporal correlation functions of single components 
of velocity differences, i.e. C{t) = {ux{t' +t)ux{t')). By 
averaging over the possible random vectors n one arrives 
at the formula 



C{t) 



1 



(16) 



i.e. there is rapid (exponential) decay with a zero-crossing 
at t* — ^TrB~^. Exponential decay and zero-crossings 
are indeed observed for the experimental data [3, H, Q 
as well as in Lagrangian DNS 13] ■ The experimentally 
observed zero of the correlation function [3, 'sf can be used 
to estimate the size of the parameter B. DNS [13, E3 
indicates that t* w 2.2r,, independent of R\, thus our 



0.95r-i. 



model parameter B is given by B •■ 

Higher-order correlation functions are of interest, too. 
For example, the correlation function of the square of 
velocity differences C{t) = {u^lf + t)'tP{t)) can be ap- 
proximated as C{t) « {l3-^{t' + t)l3-^{t')) . Clearly, by 
construction of the superstatistical model, this correla- 
tion function decays very slowly, since the process f3{t) 
evolves on a much larger time scale than the process 
u{t). This makes it clear why the correlation function 
of the modulus \u\ — \/]u]^ decays very slowly, asjt is 
indeed observed in experiments [3, |7| and in DNS [lO]. 

Finally, let us comment on the typical evolution of 
dissipation and enstrophy fluctuations in our model. A 



large value of local energy dissipation e corresponds to 
a large value of , since in our superstatistical dynam- 
ical model e ^ ^ a'^/{2j). This means the forcing 
aL{t) is strong for a while. After some short relaxation 
time of the order ~ r,,, where is the Kolmogorov 
time, this will create a large local acceleration variance 
(a^) ^ (u^)t~^ ^ P~^T~^. The term Bn x u will then 
create a lot of rotational energy (=enstrophy fl), as soon 
as \u\ has become large. Thus energy dissipation and 
enstrophy are strongly correlated in time, and enstrophy 
lags behind dissipation evolution by something of the or- 
der Tjj (the relaxation time of the system). This is actu- 
ally experimentally observed in Fig. 1 of [2^ . The peaks 
of e and fl are shifted by about half a second, which cor- 
responds to the Kolmogorov time of the system studied 
by Zeff et al. [11]. Our simple superstatistical model de- 
scribes these effects in a correct way. The above time-lag 
effect also shows up as an asymmetry of the dissipation- 
enstrophy cross-correlation function in DNS 1^, 0] ■ 

To conclude, we have demonstrated that the most im- 
portant statistical phenomena that have been experimen- 
tally reported in LT experiments so far are well repro- 
duced by a superstatistical model that can be regarded 
as a generalized Brownian motion model relevant for 3- 
d Lagrangian tracer dynamics. The model arises out of 
a physically plausible local momentum balance equation 
for the Lagrangian particle, and, compared to other re- 
cent Lagrangian models [3^ , has the advantage of being 
analytically treatable. The model naturally incorporates 
a superposition of several stochastic processes, a fast one 
for velocity differences of the tracer particle and slow 
ones for dissipation and enstrophy in the environment 
of the tracer particle. The predicted Lagrangian scaling 
exponents the 1-point and 3-point probability densi- 
ties, correlation functions, as well as the statistical depen- 
dencies between acceleration components are in excellent 
agreement with the experimental data. 
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